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We develop and analyze a theoretical model which yields the shifts and widths of Feshbach res- 
onances in an atomic waveguide. It is based on a multichannel approach for confinement-induced 
resonances (CIRs) and atomic transitions in the waveguides in the multimode regime. We replace 
in this scheme the single-channel scalar interatomic interaction by the four-channel tensorial poten- 
tial modeling resonances of broad, narrow and overlapping character according to the two-channel 
parametrization of A.D.Lange et. al. [35| . As an input the experimentally known parameters of 
Feshbach resonances in the absence of the waveguide are used. We calculate the shifts and widths 
of S-, d- and g-wave magnetic Feshbach resonances of Cs atoms emerging in harmonic waveguides 
as CIRs and resonant enhancement of the transmission at zeros of the free space scattering length. 
We have found the linear dependence of the width of the resonance on the longitudinal atomic mo- 
mentum and quadratic dependence on the waiveguide width. Our model opens novel possibilities 
for quantitative studies of the scattering processes in ultracold atomic gases in waveguides beyond 
the framework of s-wave resonant scattering. 
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PACS numbers: 32.60.-(-i,33.55.Be,32.10.Dk,33.80.Ps 



I. INTRODUCTION 



Impressive progress of the physics of ultracold quan- 
tum gases have opened new pathways for the study 
of low-dimensional few-body systems (see, for example 
P, as well as strongly correlated many-body sys- 
tems [3, 13). Specifically, it was shown that the con- 
fining geometry of atomic traps can drastically change 
the scattering properties of ultracold atoms and induce 
resonances in the collisions (confinement-induced reso- 
nances (CIRs)) (5|. The CIR for bosons has been found 
to occur in the vanishing energy collisional limit as a 
consequence of the coincidence of the binding energy of 
a diatomic molecular state with the energy spacing be- 
tween the levels of the confining (harmonic) potential 
Q. It was shown that this coincidence leads to a diver- 
gence of the effective interatomic coupling constant giD 
and to a total atom-atom reflection, which appears as a 
broad dip in the transmission coefficient T, thereby cre- 
ating a gas of impenetrable bosons @. CIRs have also 
been extensively studie d e.g . in the context of three-body 
[13, [U, and four-body |l2|| scattering in confining traps, 
fermionic p-wave scattering [l3|, distinguishable atom 
scattering [3, [l3 - IT6| or multichannel scattering [13, [3 
in atomic waiveguides. Two novel effects were predicted 
for distinguishable atoms: the so-called dual CIR yield- 



ing a complete suppression of quantum scattering jl4| . 
and the resonant molecule formation in tight waveguides 
[iGj . D-wave resonant scattering of bosons in confining 
harmonic waveguides has been analyzed very recently 
in 19]. Remarkable experimental progress has lead to 
the observation of CIRs, for identical bosons [20l - [23l | and 
fermions ^24. - 26.] . as well as to distinguishable atoms [27| . 



saeidian@iasbs.ac.irj 
melezhik@theor.jinr!ru| 



peter . schmelcher ® physnet . uni-hamburg. de | 



However, the recent experimental |23l. l26l - l28l | and theo- 
retical [13, 129I - I31I I investigations of the CIRs clearly show 
that, despite the impressive progress, the existing theo- 
retical models of CIRs need to be improved for a quan- 
titative description of the experiments in this field. The 
obvious drawback of the existing theoretical models for 
CIRs in atomic waveguides is the single-channel char- 
acter of the simple interatomic interactions employed. 
Single-channel potentials with zero-energy bound states 
were used for simulating magnetic Feshbach resonances in 
free space representing a necessary ingredient for the ap- 
pearance of the CIR in a confining trap (see, for example 
[ij, |j| ) . In the seminal work of Olshanii [3 and in subse- 
quent papers [3, HI] , including the recent works [13, [13] , 
the simple form of a pseudo-potential was used model- 
ing the interatomic interactions. More realistic potentials 
were used in our previous works M, 14 ; S[i3] as well as in 
the works of other authors [gI. Il3l. IbiL" 33l] which, however, 
all possess a single-channel character. The single-channel 
interatomic interaction approach permits to explore only 
the main attribute of the Feshbach resonances in the 3D 
free space, namely the appearance of a singularity in the 
s-wave scattering length — >■ ±cxi when the molecular 
bound state with energy Eb crosses the atom-atom scat- 
tering threshold at energy E = in the entrance channel. 
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However, other important parameters of the Feshbach 
resonance, such as the rotational and spin structure of 
the molecular bound state in the closed channel as well 
as the width A of the resonance characterizing the cou- 
pling r of the molecular state with the entrance chan- 
nel, were ignored. The main goal of the present work is 
to extend these theoretical approaches developed earlier 
for the CIRs and transverse excitations/deexcitation pro- 
cesses for collisions in harmonic waveguides Q to the case 
of a tensorial interaction as well as taking into account 
the width of the magnetic Feshbach resonances responsi- 
ble for the CIR. The parameters obtained from the exper- 
imental analysis of the magnetic Feshbach resonances in 
free space, namely the resonant energies Ec,i (or the cor- 
responding values of the field strengths Bc^i of the exter- 
nal magnetic field), the widths of the resonances Ai{Ti), 
spin characteristics and the background scattering length 
Ofjg, are used as input parameters of our model. 

It is well-known that the most adequate computational 
schemes including all the above mentioned parameters 
of the Feshbach resonances are multi-channel scattering 
P, [H, and multi-channel quantum-defect [11] ap- 
proaches providing a quantitative description for a broad 
regime of experimental parameters for the Feshbach res- 
onances in ultracold atomic gases. In the present work 
the two-channel potential scheme (35j | , permitting the ef- 
ficient modeling of the interatomic interaction near the 
known magnetic Feshbach resonances of the Cs gas in 
free space [3J], is included in our multichannel approach 
instead of the simplified interparticle single-channel 
potential we have used so far for analyzing CIRs and 
transverse excitations in atomic waveguides. With this 
approach we model s-, d- and g-wave Feshbach resonances 
in ultracold Cs (s-,d- and g- indicate here the "exit" chan- 
nels, i.e. they correspond to rotational quantum numbers 
of the molecular states in the closed channels), which 
were observed in free space scattering experiments [33 |. 
and quantitatively analyze the shifts and widths of the 
resonances in the harmonic waveguides at experimental 
conditions closed to the ones encountered in the works 
[22, 23]. We note that selected aspects of the model- 
ing of magnetic Feshbach resonances in optical traps and 
lattices were also considered in fSTI-lill. 



II. FESHBACH RESONANCE MODEL 

We consider the collision of two identical bosonic atoms 
in a harmonic waveguide. This two-body problem per- 
mits the separation of the center-of-mass and relative mo- 
tion yielding the following Hamiltonian for the relative 
atomic motion 

H{r, 9) ^ h^V^ + \^,^lp']i + V{r) (1) 

with p = rsin6' and the trap potential 1/2/xcj^p^. V{r) 
is the 4-channel interatomic potential and / is the unit 



matrix, r is the relative radial coordinate and ^ = m/2 
is the reduced mass of the atoms. 

Following the scheme suggested in [sl] for describing 
the three magnetic Feshbach resonances in an ultracold 
Cs gas let us suppose that initially the scattering atoms 
are prepared in one spin configuration |e) (called the "en- 
trance channel") and the "closed channels" |c, i) (i=l, 2, 
3) support S-, d- and g-wave molecular bound states, at 
-11. IG, 47.78G and 53.449G. The quantum state of an 
atomic pair with energy E is described as 



^ Vc,z(r)|c,i) 4^ Ve(r)|e) 



satisfying the Schrodinger's equation with the Hamilto- 
nian (1). A four-channel square- well potential 
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is employed to describe the colliding atoms in the "en- 
trance channel" |e) and the weakly-bound molecules in 
the "closed channels" |c, i) near a Feshbach resonance. 
For r <a = 'i'KT{l/4)^'^RvdWi we assume the attractive 
potential can support multiple molecular states - that is, 

VcVci 3> Eydw — ^^/2m-^^dw " ^^'^ induce Fes- 
hbach couplings between the channels [s^. The size a 
of the potential action is chosen to account for the in- 
teratomic interaction determined by the van-der-Waals 
(vdW) tail with the length Rydw = l/2(2mC6/?i^)^/^ 
and the energy scale Eydw = ^ /'^''^^dw 113 ' where Cq 



dW 

is the corresponding van der Waals coefficient and T{x) 
is the gamma function. For r > a, entrance- and closed- 
channel thresholds are set to be i? = and E = 00, 
respectively (See Fig.l in Ref. 35]). 

Such a choice of the interatomic interaction permits 
a simple parametrization of the atom-atom scattering in 
universal terms of the energy of the i-th bare bound state 
Ei, the Feshbach coupling strength of the bound molec- 
ular state with the entrance channel F^ and the back- 
ground scattering length abg, which is convenient for an 
analysis of experimental data near magnetic Feshbach 
resonances [l|, HI] . When the mixing between the closed 
channels and the entrance channel is weak and the back- 
ground scattering length \abg\ considerably exceeds the 
range of the interatomic interaction a, the s-wave scat- 
tering length a and the binding energy E'f, in free space 
are given by [ssl ]: 
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FIG. 1: (Color online) The s-wave scattering length a of \F = 
3, mp = 3) Cesium atoms as a function of the magnetic field 
B. The solid curve shows the analytical result (5) and the 
dots show the numerical result (see the text). 
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respectively. Then, assuming that the bare bound states 
can be linearly tuned magnetically by a linear Zeeman 
shift -namely, Ei — dfj,i{B — Bc^i), where 5fj,i is the rel- 
ative magnetic moment between the entrance and i-th 
channels and Bc^i is the crossing field value of the zth 
bare bound state, the scattering length a{B) can be rep- 
resented as 



abg 
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0,i 



(5) 



Here _Bo,i is the ith lowest pole of a and B* the ith low- 
est zero. The width of the ith Feshbach resonance can 
be defined as = B* — Bo,i. The binding energy of 
Feshbach molecules can be measured in the experiment 
by e.g. radio frequency and microwave spectroscopy [43l - 
|47| |. Using the fitting parameters abn, Suj, Ti and Bc^i, 
one can fit the experimental data [43l - l47| with (g]), from 
which one can calculate a{B) (3), B* and Bq ^ (5) [35*1 . 

The result of the procedure is given in the Table I 
of ref. [sll for ultracold Cs Feshbach molecules where 
abg = 1875ao and a = 95. Tag (ao is the Bohr radius) are 
also given. With these values for i?o,i and B* we show 
the scattering length a{B) as a function of the magnetic 
field B according to Eq.JS]) in Fig.l. 

The parameters atg, Sfii, Fj, -Bc,i, B* and Bo.i from [sl] 
together with a{B) defined by Eq.(5) are used for fitting 
the diagonal terms Vd and Vg in the tensor potential (2). 
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FIG. 2: (Color online) The transmission coefficient T for the 
harmonic trap with cj_l = 14.9 kHz as a function of the mag- 
netic field B for the case of the 2-channel (solid curve) and 
4-channel (dots) potential (2). 



The nondiagonal terms hfli are defined by the formulas 



r,/2 = 29'^Vc,^, tan20, = 



2hn, 



The scattering length a{B) is then calculated for different 
B and varying parameters of the potential V by solving 
the Schrodinger equation 



m (6) 



in free space {uj± — 0) with the scattering boundary con- 
ditions 

■06(1") exp{ikz} + f{k, 6) / r exp{ikr} , ipc,i{j^) 



at fcr oo for the fixed ^ (fc = ^/2fiE/h — > 0) 
psj . The diagonal matrix B in (6) is defined as Bu — 
Sfii{B — Bi) {i = 1,2,3) and Bee = 0. After separation 
of the angular part in (6) we come to the system of four 
coupled radial equations 

for the radial part 4>a{r) of the desired wave function 
IV') = EaV'a(r)la) = J2aMr)yio,o{r)\a), where a = 
{e,i = 1,2,3}, {le = OJi = 0,^2 = 2,^3 = 4) and the 
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matrix elements Vai3{r) are defined by Eq.(2). The cen- 
trifugal barrier — ^ '2^°^"'""'^'* C'') niodels at r — )■ the 
correct asymptotic behavior 4>i{r) ^ r''"*"^ of the molecu- 
lar bound states |c, i) in the closed channels which couple 
to the entrance s-wave (Zg = 0) channel |e) by the nondi- 
agonal terms V^0(r)(a 7^ (3). 

By varying the 14,^, T4 and $7^ we obtain an excellent 
agreement of the calculated s-wave scattering length a{B) 
with the analytical results from [35] for Cesium atoms in 
the hyperfine state \F — 3,mp — 3) for the considered 
magnetic field regime — 40G < B < 60G (see dots in 
Fig.l). In this regime, we observe three resonance terms, 
which correspond to the coupling to the s-, d-, and g-wave 
molecular states. 

Next we analyze the scattering properties of the s-, d- 
and g-wave magnetic Feshbach resonances in harmonic 
waveguides by integrating the Schrodinger equation (6) 
for ujj_ 7^ with the scattering boundary conditions 

-0e(r) = (cos(fco2:) + /eexp{jfco I Z |}) $o(/5) , ^Jc.iir) 

(8) 

at I z I = I r cos | — 00 adopted for a confining trap 
0. fe{E) is the scattering amplitude, correspond- 
ing to the symmetry with respect to the exchange 
z ^ —z (we consider collisions of identical bosonic 
Cs atoms), $o(p) is the wave function of the ground- 
state of the two-dimensional harmonic oscillator and 
ko = ^2[i(E — hjj±^/h = ^J2^E\\/h. In the presence 
of of the harmonic trap (a;j_ 7^ 0) the problem (6), (8) 
becomes non-separable in the plane {/?, z}, i.e. the az- 
imuthal angular part is separated in the wave function 
IV') and Eq.(6) is reduced to the coupled system of four 
2D Schrodinger-type equations. To integrate this cou- 
pled channel 2D scattering problem in the plane {r, 6*} 
we have extended the computational scheme 3- 

The computations have been performed in a range of 
variation of uj±_ close to the experimental values of the 
trap frequencies ~ 27r x 14.5kIIz [2^. We have integrated 
Eq.(6) for varying B and fixed longitudinal colliding en- 
ergy i?| I ^ E—fiLOj_. In the main part of computations the 

energy was chosen very low — 1.0 x 10~^^(^^) 

^. (fco = 1.0 X 10"^(^) 0) to have a possibility 
for direct comparison with existing pseudopotential es- 
timates obtained in zero-energy limit. The integration 
was performed in the units of the problem leading to the 
scale transformation: r— >■ =,_E— > -S-, V ~> and 

uj^-^ with £^0 = and = ^. 

III. RESULTS AND DISCUSSION 

A. Transmission coefficient 

In Fig. 2 we present the transmission coefficient T{B) = 
|H-/e(i3)P as a function of B calculated for the harmonic 
trap Lo±_ = 14.9kHz with the two-channel (i=l) and four- 
channel (i=l,2,3) character of the tensorial interaction 
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FIG. 3: (Color online) The coupling constant gm as a func- 
tion of the scattering length a in free space in the region near 
the s-wave Feshbach resonance, calculated for the harmonic 
trap with lj± — 14.9kHz. Dots indicate the numerical data 
and the solid curve the analytical results with the formulas 
derived in the pseudopotential approach jsj . 

V{r) (2). The two-channel potential V{i = 1) supports 
only one broad s-wave Feshbach resonance while the four- 
channel potential V{i = 1,2,3) supports all three s-, d- 
and g-wave magnetic Feshbach resonances in free space 
(a;_L = 0) at the fields Bo,i = -U.IG, Bo,2 = 47.78G and 
-So, 3 = 53.449G, respectively All these resonances 

are also developed in the calculated curve T{B) of the 
transmission coefficient in the trap ^ 0). 
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FIG. 4: (Color online) The transmission coefficient T{B) as a 
function of the strength of the magnetic field for several trap 
frequencies uj±. Here the s-wave scattering length a{B) from 
Fig.l is also given (solid line). 
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FIG. 5: (Color online) The dependence of the s-wave scat- 
tering length a{B) on the width a± of the waveguide at the 
points Bmin of the minimum of the transmission coefficient 
T{a{B),a±) (see Fig. 4). The dots, pluses and cycles cor- 
respond to the calculated points near the s-,d- and g-wave 
magnetic Feshbach resonances, respectively. The solid curve 
corresponds to the formula a = a±/C 0]. 

First, we have analyzed the region of -B ~ 
— 11.1 G near the s-wave resonance. In Fig. 3 we 
show the calculated ID coupling constant giD = 
linifcn_>o Re{fe{ko)}/ Im{fe{ko)}ko/ fj. [6| as a function of 
the scattering length a{B) in this region. By tuning 
B in the interval -60G < B < 40G , the s-wave 
scattering length a{B) changes from —oo to +oo. At 
B = — 49.78G corresponding to the point a(B) = a±/C 
(where C = 1.4603.. and aj_ — y/h/^ujj_) of the appear- 
ance of the CIR [5], the coupling constant gioio-/ a-i) di- 
verges and the behavior is in very good agreement with 
our previous computation of (?iD(a/a^) performed with 
the single channel screened Coulomb potential [Sjj and 
the formulas giD = 2fta/(uaj_)/(l — Ca/a±) derived in a 
pseudopotential approach [5||^ 

Let us next analyze the region of the d-wave magnetic 
Feshbach resonance near the point 47.78G (see Fig.l) - 
the region of major experimental interest due to atomic 
loss and the formation of Cs molecules in the confined ul- 
tracold gas of Cs atoms [13, HH . The results of our com- 
putations are illustrated by the curves T{B) calculated 
for different uj± (Fig. 4) corresponding to the transverse 
frequencies of the optical trap being used in the exper- 
iment [23| . It is clearly shown that the position Bmin 
of the transmission coefficient T{B) minimum is depen- 
dent on the trap width a± = yjhj {^uj±) and the corre- 
sponding scattering length a{B„iin) at the point Bmiti is 
accurately described by the formulas a{Bmin) = cl^/G 
obtained by Olshanii ^ for the position of the CIR in 
a harmonic waveguide. The latter fact is illustrated in 
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FIG. 6: (Color online) (a) The dependence of the width F* 
on the longitudinal momentum fco near the d-wave Feshbach 
resonance. The solid line has been obtained via Eq.(lO), solid 
circles indicate the widths r*(fco) extracted from the numer- 
ically calculated r(_B,fco).(b) The transmission coefficient T 
as a function of the magnetic field B calculated for a few ko. 
For both subfigures iJx = 14.9kHz. 



Fig. 5. Here the results of the calculation of the depen- 
dence of a(Bmin) on a± for s-wave and g-wave Feshbach 
resonances are also given. This analysis clearly demon- 
strates that the law a = a±/C for the position of the 
CIR in a harmonic waveguide is fulfilled with high ac- 
curacy for the Feshbach resonances of different tenso- 
rial structure, although the law was initially obtained at 
zero-energy limit in the framework of an s-wave single- 
channel pseudopotential approach to the interatomic in- 
teraction |5|. We also see that the positions B^ax of the 
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FIG. 7: (Color online) The dependence of the width F* on 
the trap frequency tJx near the d-wave Feshbach resonance at 
fco=2.0899x 10~* (^). The solid curve has been obtained via 
Eq.(lO), solid circles indicate the widths r*(a)x) extracted 
from the numerically calculated T{B,ijj±). The analytical 
(dashed line) and numerical (squares) results for the right 

CIR 



half width l/2Fg}^ versus bj± are also given. 



maximums of the coefficients T{B) calculated for differ- 
ent uj± are independent of uj±_ and coincide at the point 
B2 = 47.944G of the zero of the scattering length in free 
space {a{B2) = 0). This fact is in agreement with the 
analytic result T —\ 1 + fe p— 1 obtained in the pseu- 
dopotential approach at a — > 0. The same behavior of the 
T{B) coefficients has been found in the vicinity of 18. IG 
and 53.46G (points B^ and B^ of s-wave scattering length 
zeros near the s- and g-wave Feshbach resonances). 



B. Width of the resonant enhancement of the 
transmission 

While the coefficients T{B) in the vicinity of the points 
B* (positions of the zeros of the scattering length a{B)) 
show a resonant behavior one can define the width F* 
of this resonance as the width at half T maximum. By 
using the formulas for the even scattering amplitude in 
the trap [11, Ei] 
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ID 



valid near the CIR in the zero-energy limit and the def- 

r3 (B-Bl) 



mitions a ID 
we obtain 



^(C-^) anda = a,,n-=i(B3Bi 
4af,g7ia^fcoAi 



(9) 



where 7^ = Ilj^^i (b'-Bq \ • '^^^ above formula is valid 
for the condition F* ^ B* of narrow resonances, which 



is fulfilled with high accuracy for the d- and g-wave res- 
onances and less accurately for s-wave resonance (see 
Fig. 2). In the low energy limit fco — > the expression 
is reduced to 



r* = A. 



' abgli 



= A. 



(10) 



where A^ — B* — Bg^i is the width of the Feshbach reso- 
nance in free space and the dimensionless 7^- factor w i 
for the case of d- and g-wave resonances and ~ 1 for s- 
wave resonance. Figs. 6 and 7 demonstrate the very good 
agreement of these formulas with the numerical computa- 
tions. Actually, in Fig. 6(a) one can see in a broad range 
of ko variation the very good agreement of the width 
F*(fco), extracted from numerically calculated T{B,ko), 
with the linear functional dependence on the longitudi- 
nal momentum kg, following from Eq.(lO). Also, the 
F*(a;x) has inverse dependence on ujj_. It is confirmed 
by results of numerical calculations of F*(a;_L) given in 
Fig. 7. Thus, Eq.(lO) can be used for extracting impor- 
tant information about the analyzed system. Indeed, by 
measuring the width F* one can extract from Eq.(lO) 
the longitudinal momentum ko = ^J2^E\\/h (longitudi- 
nal colliding energy E\\ ) and estimate the "longitudinal" 
temperature of the atomic cloud in the trap. This ex- 
pression also shows that one can control the width F* 
of the resonance by varying the trap frequency In- 
creasing of cjj_ leads to a narrowing of the resonance (see 
Figs. 4 and 7), the effect of which can be used experimen- 
tally. Fig. 6(b) demonstrates the stability of the position 
of the CIR (the minimum of the transmission coefficient 
T) with respect to fco variation and the linear growing 
of the transmission coefficient T with fco increasing. It 
is shown that already at very low fco ^ 10"^ (^) the T- 
coefficient becomes large enough and consequently can 
be experimentally "visible" near CIR. 

To conclude this subsection we also give the expression 
for the position of the CIR with respect to the resonance 
at a = (to the point B*) 



^CIR 



1 - C7,^ 



(11) 



Since the resonance maximum at B* corresponds to the 
unit value T{B*) — )■ 1 and the minimum of the CIR ap- 
proaches to zero T{B* — A^^^) 0, we can define the 

"right" half width of the CIR as the distance of 

the zero of the coefficient T to the right situated point 
where T = i. Using the formulas for A^^^^ and Eqs.(lO) 
for F* we obtain 



1. 



CIR 



^CIR 



1 



a^. fco 



(fco->0) 



'1 



C7.ff 2a,g7,J ■ 
(12) 

The above expression describes the narrowing of the CIR 
width with decreasing trap frequency opposite 

to F* (see Fig. 4). This expression might also be of inter- 
est to the experimental analysis of the CIRs. 
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FIG. 8: (Color online) The relative populations Pi/Pe (13) of 
the molecular states i — s,d and g calculated as a function 
of B near the d-wave Feshbach resonance in Cs for the pair 
collisions in free space (a) and in the harmonic waveguide with 
ojx = 59.6kHz (b). The s-wave scattering length a{B) in free 
space from Fig.l is also given. 



C. Population of molecular states in free space and 
in the waveguide 

We have also calculated the relative populations Pi /Pg 
of the molecular states |c, i) in the process of pair atomic 
collisions for harmonic traps of different frequencies uj± 
as well as in free space — 0). The populations Pi are 
defined as 



\^Jc^,{r,9)\^r^drsm0d0 



(13) 



where integration over the infinite region < ^ < oo 
gives a convergent results due to the decaying tails of the 
molecular bound-state wave-functions "00,4 (i") — >■ at dis- 
tances of the order ~ a in the closed channels |c, i). The 
population Pg is defined in a region near the origin r — >■ 
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FIG. 9: (Color online) Schematic picture (presented in [a, [131 
for a single-channel potential model of interatomic interac- 
tion) of the bound state E\, of two atoms as a function of 
ax I a in free space (solid line) and in the harmonic waveguide 
(dashed line). The first excited resonance state Er is also 
presented. 



of the molecular dimension ~ a in the entrance channel 
by using Eq.(13) where the upper limit of the integration 
over r is a. The relative populations Pi / Pe of the molec- 
ular states for ultracold atomic collisions in free space 
and in the harmonic waveguide with a;_L = 59.6kHz are 
shown in Fig. 8 in the region of B near the d-wave Fesh- 
bach resonance. In free space the pair collision leads to 
a relatively low population of the d-wave molecular state 
P2/Pe ~ 10~^ — 10~^. The populations of non-resonant 
s- and g-wave states are Pi ^ P2 and P3 <C P2 essen- 
tially suppressed here. The dependence of the d-wave 
molecular state population on B repeats the dependence 
on B of the population of the region ^ a in the entrance 
channel P2(P)/Pe(P) ^ const. However, in the confined 
geometry of the waveguide (cjj_ 7^ 0) the pair collision 
leads to a resonant enhancement of the relative popu- 
lation P2(P)/Pe(P) of the d-wave molecular state near 
the point P2 = 47.9446 where the free space scattering 
length is zero a(i?2) = 0. In the waveguide the point of 
appearance of the bound state is shifted from the posi- 
tion defined by l/|a| — >■ for free space scattering (see 
point A in the illustrative scheme of the bound and reso- 
nant states of an atomic pair given in Fig. 9) to the point 
l/|a| — 00 (the point B in Fig. 9). This is why we observe 
in Fig. 8(b) the strong resonant enhancement of the pop- 
ulation P2(P2) at the magnetic field B^ = 47.944G. The 
populations Pi/Pg and P^/Pe of other molecular states 
also show some enhancement with respect to the entrance 
channel due to the coupling of the states with the "reso- 
nance" channel |c, 2), which, however, are a few orders of 
magnitude less than the enhancement of the population 
PijPe- We do not observe a resonant behavior of P2(P) 
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at the point B = 47.57G of the CIR (see Fig.4 and the 
point C in Fig. 9). We suspect that this is the case due 
to the rather weak coupling in our model potential V{2) 
between the entrance channel |e) and the resonant state 
in the closed channel |c, 2) for stimulating considerable 
transition to the molecular state |c, 2) in the closed ex- 
cited channel of the waveguide. 

In Figs. 10 and 11 we present the results of our com- 
putation of the molecular state populations near the s- 
and g-wave Feshbach resonances -11. IG and 53.449G, 
respectively. We observe the qualitatively analogous 
effect of resonant enhancement of the relative popu- 
lations Pi{B)/PeiB) and P3{B)/Pe{B) at the points 
B* = 18. IG and ^3 = 53.457G in the harmonic waveg- 
uide, which correspond to the positions of the zero of 
the s-wave scattering length in the vicinity of s- and 
g-wave Feshbach resonances, respectively. Here we no- 
tice Pi{B)/Pe{B) ~ const, i.e. the same dependence 
on B of the s-wave molecular state population as the 
population of the region < a in the entrance channel 
Pi(B)/Pe(B) ~ const (see Fig.lO(a)). In Fig.ll(a) we 
observe the resonant enhancement of the relative popu- 
lation P2,(B)/Pe{B) of the g-wave molecular state near 
the Feshbach resonance in free space at i?o,3 = 53.449G 
where l/|a| — 0. A similar resonance enhancement 
near the point l/|a| — >■ of the appearance of the 
near threshold resonance or weakly-bound state in free 
space (see point A in Fig. 9) is observable via the values 
Pi{B),P2{B) and Pe{B) but, because Pi{B)/Pe{B) ~ 
const and P2{B)/Pf.{B) const, we do not observe 
this effect in the relative populations Pi{B) / Pe{B) and 
P2{B)/PeiB) in Figs.lO(a) and 8(a). 



IV. CONCLUSION 

We have developed a theoretical model for a quan- 
titative analysis of the Feshbach resonance shift and 
width induced by an atomic waveguide. It is based 
on our multichannel approach for confinement-induced 
resonances and atomic transitions in the waveguides 
in the multi-mode regime ^J. In this scheme the 
single-channel (scalar) interatomic interaction is replaced 
by a four-channel (tensorial) potential modeling reso- 
nances of different structure according to the two-channel 
parametrization of A.D.Lange et. al. [s^. The experi- 
mentally known parameters of the Feshbach resonance in 
the absence of the waveguide are used as an input in our 
approach. We have calculated the shifts and widths of s-, 
d- and g-wave magnetic Feshbach resonances of Cs atoms 
emerging in harmonic waveguides as CIRs and resonant 
enhancement of the transmission at zeros of the free space 
scattering length. 

In particular we find that the relationship a — a±/C 
for the position of the CIR in a harmonic waveguide is 
fulfilled with high accuracy for the Feshbach resonances 
of different tensorial structure which holds in spite of 
the fact that this property was originally obtained in the 
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FIG. 10: (Color online) The relative populations Pi/Pe (13) 
of the molecular states i = s,d and g calculated as functions 
of B near the s-wave Feshbach resonance for Cs for the pair 
collisions in free space (a) and in the harmonic waveguide with 
Lj± = 59.6kHz (b). The s-wave scattering length a(B) from 
Fig.l is also provided. 



framework of a s-wave single- channel pseudopotential ap- 
proach Note, that this property was experimentally 
confirmed for d-wave Feshbach resonances in a gas of 
Cs atoms [23j . The maximum of the transmission, cor- 
responding to the zero of the scattering amplitude, is 
shown to be independent of the trap field strength and 
then again corresponds to the zero B* of the s-wave scat- 
tering length a{B) in free space. In a nutshell, the Fes- 
hbach resonance in free space develops in the harmonic 
waveguide into a minimum of T (position of CIR), de- 
fined by the formulas a = a±/C, and a maximum, coin- 
ciding with the position of zero of the s-wave scattering 
length a. The 'distance' between these extrema is equal 
to A^^^^ = Ai/(1 - Cjiabg/a±_). 

We have derived expressions for the widths F*, ^F^^^ 
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side half-width of the minimum of the T-coefScient i.e. at 
the position of CIR and confirmed its vahdity by numer- 
ical results for fco ~^ 0. By measuring the width T* one 
can, in principle, extract from these expressions the longi- 
tudinal collision energy and estimate the "longitudinal" 
temperature of the atomic cloud in the trap. In other 
words the width of the atomic loss resonance observed 
in the experiment [23| at the point of CIR might contain 
important information about the longitudinal atomic mo- 
mentum ko and the temperature of the gas. It also shows 
that one can control the width T* of the resonance at 
a{B*) = by varying the trap frequency uj±. An in- 
crease of uj± leads to a narrowing of the resonance, an 
effect which could potentially be used experimentally. 

Finally, the molecule formation rates in a waveguide 
show an enhancement for the case of a corresponding 
zero of the s-wave scattering length a{B*) — 0. We have 
shown that the positions of these resonances are stable 
w.r.t. the variation of the confining frequency a;_L of the 
waveguide. 

Our model adds to the possible studies of scattering 
processes of ultracold atomic gases in waveguides beyond 
the framework of s-wave resonant scattering. Our model 
might be extended to the cases of fermions or distin- 
guishable atom scattering, including transverse excita- 
tion/diexcitation processes [s*]. It permits also for the 
investigation of other trap geometries [48] and more re- 
alistic interatomic interactions. 
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